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A unified response theory for the time-resolved nonlinear light generation and two-photon pho- 
toemission (2PPE) from metal surfaces is presented. The theory allows to describe the dependence 
of the nonlinear optical response and the photoelectron yield, respectively, on the time dependence 
of the exciting light field. Quantum-mechanical interference effects affect the results significantly. 
Contributions to 2PPE due to the optical nonlinearity of the surface region are derived and shown 
to be relevant close to a plasmon resonance. The interplay between pulse shape, relaxation times 
of excited electrons, and band structure is analyzed directly in the time domain. While our theory 
works for arbitrary pulse shapes, we mainly focus on the case of two pulses of the same mean fre- 
quency. Difficulties in extracting relaxation rates from pump-probe experiments are discussed, for 
example due to the effect of detuning of intermediate states on the interference. The theory also 
allows to determine the range of validity of the optical Bloch equations and of semiclassical rate 
C ' equations, respectively. Finally, we discuss how collective plasma excitations affect the nonlinear 

^ optical response and 2PPE. 
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■ I. INTRODUCTION 

o\ 

, During the last decade time-resolved spectroscopy of condensed-matter systems has become a very active area 

■ of experimental research [1-22]. This is mainly due to the progress in experimental technique, in particular the 
T— I ' ability to create ultra-short laser pulses with a duration of the order of a few femtoseconds [23]. Since this is 

similar to the relaxation times of excited electrons and collective excitations in solids, these experiments allow to 
study non-equilibrium physics, e.g., the time evolution of excited electrons before and during thermalization. Of 
particular interest are non-linear techniques such as time-resolved sum-frequency generation (SFG) and two-photon 
photoemission (2PPE), which are sensitive to excited electron states [24]. A theoretical understanding of these 
processes is crucial. Petek and Ogawa [21] noted in 1997 that a theory for time-resolved 2PPE is still lacking, and, 
^ ■ despite the efforts of many theorists, much remains to be done. The situation for SFG is similar. The construction 
Q ' of such a theory is a formidable task — the main problems are (a) the desription of the time-dependent response and 
O (b) the treatment of the surface. Our main concern is with the first point. A simplified description of the surface 
^ _ using Fresnel factors has been employed successfully to describe SFG from metals [25-28]. A detailed discussion of 
boundary conditions at the surface, focusing on the nonlinear optical response of magnetic systems, can be found in 
Ref. [29]. 

In the present paper we discuss the electronic processes taking place during time-resolved SFG (in particular 
second-harmonic generation, SHG) and 2PPE and derive the dependence of the SFG light intensity and the 2PPE 
photoelectron yield on the time dependence of the exciting laser field. We show that most effects observed for time- 
resolved 2PPE appear similarly for SFG, such as their dependence on energy relaxation, dephasing, and detuning 
of intermediate states. Other examples are the enhancement of the response due to collective excitations and the 
sensitivity regarding the ultra-fast spin-dependent relaxation. We develop a unified time-dependent response theory 
for SFG and 2PPE, starting from the self-consistent field approach of Ehrenreich and Cohen [30, 31], which can be 
applied to specific materials described by their band structure, relaxation rates, and dipole matrix elements. For 
illustration, we apply the theory to a generic tight-binding model for a metal to study interference effects in both 
pump-probe single-color SFG and 2PPE and their dependence on relaxation rates and detuning. We exhibit the 
strong similarities between both methods. 

In SFG [1-5] electrons are excited by absorbing two photons and they subsequently emit a single photon at the sum 
frequency. In Fig. 1 we illustrate the type of process yielding SFG. For simplicity we talk about SFG in the following, 
although dij(ference-frequency generation is automatically included in our theory. Time-resolved measurements [1-5] 
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FIG. 1: Simplified representation of sum-frequency generation (SFG). Ef is the Fermi energy and -Bvac is the vacuum energy. 
In SFG two photons of frequencies cji and cj-z are absorbed by electrons in states |1) and |2) and a single photon of frequency 
wi + ui2 is emitted due to the electronic transition |3) The two photons may be provided by one or two laser pulses. 

Note, whether the two photons are predominantly absorbed at nearly the same time or with some delay depends on the shape 
and width of the pulse(s). A more careful analysis of the absorption process on the basis of response theory is given in Sec. II A. 



usually employ the pump-probe technique, where two laser pulses of the same (single-color) or different (two-color) 
frequency are applied with a time delay AT between them. This time delay controls the time between the two 
absorptions and thus the relaxation dynamics of the electron in the intermediate state |2) is crucial, see Fig. 1. SFG 
is strongly surface-sensitive, since the SFG response of the bulk of an inversion-symmetric crystal vanishes in the 
dipolc approximation. The inversion symmetry can also be broken by nanostructures. The most important case of 
SFG is second-harmonic generation (SHG), where the electrons are excited by approximately monochromatic light 
of frequency lo and light of frequency 2u} is detected. Note, in the case of ultra-short laser pulses the spectrum is 
necessarily broadened and a full treatment of SFG is required even for these single-color experiments. Also note that 
a single laser pulse, depending on its duration and shape, involves time-delayed absorptions. 

Time-resolved 2PPE experiments of metal surfaces [6 18] as well as of clusters [19, 20] employing the pump-probe 
technique have been performed more often than time-resolved SFG. Reviews can be found in Rcfs. [21] and [22]. 
Figure 2 shows a sketch of the processes yielding 2PPE. An electron is excited above the vacuum energy E'vac due 
to the absorption of two photons. The interplay between the relaxation of the electrons in intermediate states and 
the time between the two absorptions will determine the resulting photoelectron current. The probability of electrons 
above the vacuum level actually leaving the solid is also crucial. The limited mean free path of the electrons makes 
photoemission surface-sensitive, but in general less than in the case of SFG. In both SFG and 2PPE interference 
effects [1, 3, 8, 11-13] appear, which our theory allows to study. Of course, these interference effects are expected to 
depend on the pulse shapes. 




FIG. 2: Simplified representation of two-photon photoemission (2PPE), where Ep is the Fermi energy and -Bvac the vacuum 
energy. Here, two photons of frequencies uji and L02 (out of the same or different pulses) are absorbed by electrons in states |1) 
and 1 2), respectively. Compared to SFG, Fig. 1(a), the excitation energy is now so large that electrons are excited above Svac 
and can leave the solid. The open arrow denotes electrons leaving the crystal. Note, SFG is also possible due to a transition 
|3) — » |1). However, the SFG intensity may be small, since it involves more dipole matrix elements, as we discuss below. 



The response theory presented here goes beyond previous theoretical treatments of ultra-fast processes [32] in 
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SFG and 2PPE in metals, which mainly fall into four classes: (a) density functional theory and approaches based 
thereon [28, 33 38], (b) rate equations [15, 39, 40], (c) optical Bloch equations [9, 12, 43], and (d) perturbative methods 
[26, 31, 44, 45]. At first, density fmictional theory has been applied in the time-dependent local-density approximation 
for jellium models [33-36] . In the jellium approximation one ignores the potential of the ion cores and, consequently, 
any band-structure effects. Thus this approach is not suitable if single bands or surface states or quantum-well states 
in thin films are important. On the other hand, collective excitations are usually described rather well [46]. Going 
beyond the jellium model, Luce and Bennemann have employed the local density approximation to calculate dipole 
matrix elements as they enter also in our approach [28]. Additionally taking excited states into account within the 
GW approximation, Schone et al. [37] have calculated electronic lifetimes. Hole dynamics have also been studied 
with density functional methods [38]. 

However, one would like to gain more general physical insight than the numerical results can provide. To this end 
one may consider rate equations for the occupation of excited states, e.g., the Boltzmann equation [15, 39, 40]. This 
approach allows to incorporate important effects such as secondary electrons due to relaxation from higher-energy 
states and to Auger processes as well as transport into the bulk [15, 39, 40]. However, rate equations neglect the 
electric polarization of the electron gas, its dephasing, and any quantum-mechanical interference effects, resulting 
from the superposition of the laser field and the induced fields. To include these effects one has to solve the equation 
of motion for the entire density matrix p, not only for its diagonal components, i.e., the occupations. This can be done 
in response theory. Its simplest form yields the optical Bloch equations: The system is modelled by a small number 
of levels and the von Neumann equation of motion (master equation) for the density matrix is integrated numerically 
[9, 12, 43]. However, this approach is limited to a small number of levels so that a realistic band structure cannot be 
described. Furthermore, many-particle effects like collective excitations are not included. 

On the other hand, the response theory presented here does include the band structure and collective excitations. It 
generalizes the theory of Hiibner and Bennemann [31] to SFG due to incident light of arbitrary time dependence and 
spectrum. The previous theory [31] has been used successfully for SHG from metal surfaces, thin films, quantum wells, 
and metallic monolayers due to continuous-wave, monochromatic light [26-28, 31, 47-50]. However, the dependence 
of SHG on the pulse shape and the effect of energy relaxation and dephasing were not discussed. We also derive the 
response expressions for time-dependent 2PPE within the same framework. Since our theory is explicitly formulated 
for continuous bands, it can also serve as a basis for the discussion of the averaging effects due to bands of finite 
width discussed in a more heuristic framework using optical Bloch equations for discrete levels in Ref. [51]. Since the 
full time or frequency dependence is included, effects of frequency broadening of short pulses and of finite frequency 
resolution of the detector (for SFG) [51] are easily studied. 

Our theory employs a generalized self-consistent-field approach [30, 31], which is equivalent to the random-phase 
approximation (RPA) [52-56] . We employ the electric-dipole approximation, which is valid for small wave vector q of 
the electromagnetic field and has been used successfully to describe SHG from metal surfaces [26 28, 47 50, 57]. This 
is reasonable, since the skin depth, which is the length scale of field changes, is about one order of magnitude larger 
than the lattice constant. One has to take care in interpreting SFG experiments for inversion-symmetric crystals, 
since the surface contribution only dominates over higher multipole bulk contributions for surfaces of low symmetry 
[26, 31]. Similar in spirit to our response theory, Ueba [58] has studied continuous-wave 2PPE from metal surfaces, 
Pedersen et al. [44] have considered continous-wave SHG from metallic quantum wells, and Shahbazyan and Perakis 
[45] have developed a time-dependent, but linear response theory for metallic nanoparticles. 

It is important to understand that at the surface of a metal, in thin films, and in nanostructures the light couples to 
collective plasma excitations. The field within the metal is of course not purely transverse [41, 42]. Its transverse and 
longitudinal components couple with the conduction electrons to form plasmon-polaritons and plasmons, respectively 
[42] . The (longitudinal) plasmon modes only decouple from the applied field for a structureless jellium model of the 
solid [41, 42]. However, we consider a more realistic model that incorporates the crystal structure. Also, we will see 
that the induced nonlinear polarization couples to (longitudinal) plasmon modes. 

On general grounds one may expect that the discussion of the intimate relationship between 2PPE and SFG also 
helps to understand the dependence of 2PPE on light polarization. It has been shown that the light-polarization 
dependence of SFG is important for the analysis of the electronic structure and magnetism [59] . 

The organization of the remainder of this paper is as follows: We first summarize the response theory for SFG and 
2PPE in Sec. II. This lays the ground for our discussion in Sec. HI of time-dependent SFG and 2PPE. Details of the 
response theory are given in appendices A and B. 
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II. RESPONSE THEORY 
A. Sum-frequency generation 

We first outline the response tlieory for SFG. We consider a semi-infinite solid with single-particle states |k||Z) with 
energies iJk|,; described by the momentum k|| parallel to the surface, which is assumed to be perpendicular to the 
z direction, and a set of additional quantum numbers I. For bulk states, which may be affected by the surface but 
are not localized close to it, the composite band index is I = (fc^,;^, a), where is the z momentum component in 
the bulk, v is a. band index, and a is the spin quantum number, kz has a continuous spectrum. On the other hand, 
for states localized at the surface, I is discrete. Examples are image-potential states, adsorbate states, quantum-well 
states in a thin overlayer. and proper surface states. 

Part of the electron-electron interaction is included by the self-consistent- field approximation or RPA [30, 52]. The 
remaining electron-electron scattering is approximately taken into account by inserting phenomenological relaxation 
rates [60] into the single-electron Green functions and by shifting the band energies Ey^^i [61]. We assume that i^ky; 
are quasiparticle energies containing these shifts. Note, the clectron-phonon interaction only becomes relevant on 
longer time scales and is not considered here [32]. Also, miraband contributions to the response are not considered 
for simplicity, which is reasonable at optical frequencies. 

The electrons are coupled to the effective electric field E within the solid through a dipolar interaction term (for 
simplicity we assume that the dipole coupling dominates). The optically induced polarization P within the solid is 
expanded in orders of the electric field E. The linear response is given by 



(1) 



where Xij is the linear susceptibility, q = {(i\\,qz), and qy = qy due to conservation of momentum parallel to the 
surface. Summation over repeated indices is always implied. T?he non-conservation of is explicitly taken into 
account. 

We assume throughout that the photon momentum q is small compared to the dimensions of the Brillouin zone 
and that the band energies, relaxation rates, and transition matrix elements change slowly with momentum so that 
the difference between the parallel crystal momentum of an electron before and after the interaction, ky and k'| , 
respectively, can be ignored. If we further neglect the frequency dependence of the transition matrix elements the 
self-consistent-field approach gives the time-dependent linear susceptibility 

Xy(q,q';t-t') = ^^e(t-i') 

^ E E ^k||/';k||;(-9.) ^'k||/;k,|!' Wz) 
k|, W 



X [f{Ek..i') - fiEk::i)] exp 



n 

xexp[-rk||;;k||i'(i-i')], (2) 

where v is the volume of the system. Note, the last two factors explicitly describe the oscillations and decay of the 
linear induced polarization. In the dipole approximation the transition matrix elements are 

Dk||;;k||H9.) = (k||«|r|k||0- (3) 

The matrix elements are given without approximations in App. A. The linear susceptibility is represented by the 
usual Feynman diagram shown in Fig. 3. 

The finite lifetime of electrons due to their interaction enters Eq. (2) through the dephasing rates rk||;;kj|r, which 

describe the decay of the superposition of states |k||Z) and |k||Z') and thus of the polarization. The change of the 

occupation of states is described by the energy relaxation rates rk||i;k||i = ''"i^i^;; where ti^^^i are the lifetimes. rk||;;k||; is 
the rate of spontaneous transitions out of the state \k\\l). Since the depopulation of the states or |kj|i') certainly 
leads to the destruction of the polarization, the dephasing rates can be expressed in terms of the lifetimes as [61] 

J-k||i;k i'— „ "I" k||J;k;,i" 
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FIG. 3: The Feynman diagram for the hnear susceptibihty x relating the linear polarization P'^-* to the cfTective electric field E, 
cf. Eq. (1). The solid lines in the diagrams are to be understood as electronic Matsubara-Green functions containing rela^Kation 
rates T. The dots (•) denote dipole matrix elements D. 

where T^^ describes additional dcphasing. 

The induced second-order polarization is given by 

PF\t) = j dh dh xfkit - *i - *2) Ej{U) Ek{t2). (5) 

The second-order susceptibility x^^"* depends only on two time differences due to homogeneity in time. Obviously, 
\ti — is the time interval between the two absorptions. For a single laser pulse this interval is controlled by the 
pulse width. For two pulses we expect a contribution for \ti — ^2! of the order of the pump-probe delay time AT. 
Note, the light polarization is characterized by the components Ej . 

To express the electric field E within the solid in terms of the applied external light field Eias and similarly the 
electric field Eout of the outgoing light in terms of the polarization P one should employ Frcsncl formulae, which are 
also of importance for the coupling of the light to collective excitations, as we discuss below. We do not present the 
Presnel formulae here, since they can be found in the literature [25, 26]. See also Refs. [42, 76] for effective Presnel 
factors for systems of several layers, such as the important case of a coupling prism separated from the metal by a 
thin layer of air or vacuum [77]. Of course, it would be of interest to repeat Frcsnel's analysis for SFG, in particular 
deducing phase shifts etc. 

Equation (5) is the basis for time-dependent SFG. Clearly the pulse shape of the applied light described by E(i) 
affects the induced polarization P'^^ (t). Note, for simple pulse shapes (Gaussian, Lorentzian, rectangular) it is possible 
to evaluate the integrals in Eq. (5) further. The light polarization dependence is controlled by the symmetries of the 
tensor xijl- T^^*^ symmetries of x*^^' for magnetic and nonmagnetic crystals under monochromatic light have been 
discussed in Ref. [78]. They are determined by the symmetry operations that leave the particular surface invariant. 
These symmetry arguments are unchanged for general time dependence of the applied laser field. 

The intensity of SFG Hght is I^'^\t) oc [-Eoui(^)]^ °^ [P'^^)(f)]^. So far, typical experiments do not resolve the time 
dependence of the intensity, but measure the time-integrated SFG yield 



J(2) = j dtl^^\t) (X J dt [S^^^it)^'' cx j dt [P(^H*)]'- (6) 



For simplicity we here sum over polarization directions. 

Time- resolved SFG may be performed by measuring l'^'^\AT) as a function of the time delay AT between the 
applied field pulses. Omitting surface effects (Fresnel factors) to emphasize the structure of the results, the yield can 
be written as 

J(2) (AT) Gc j dt dtx dt2 dh dU x\fk{t -ti,ti-t2) 

X XuLit - *3, h - U) Ejih) Ek{t2) Eiits) EmiU), (7) 

which is of fourth order in the incoming light field and thus of second order in its intensity. As mentioned above, the 
typical time differences dominating the response are controlled by the delay AT, besides the pulse durations. 
It is useful to write the second-order polarization P'^^' also in frequency space, 

PPiu;) = J d^'xfui.'^,J)Ei{J)Eu{uj-J), (8) 
where xgi(i - hM - ^2) = / dw' e-^"(*-*^) g-zc^'(t2-ti) ^f^^{uj^J) or 

xgi('^,'^') = 4^/ dtdi'e^-*e-*'xgi(t + i',-t')- (9) 
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Note that we employ the convention of Eq. (8) in Ref. [31] for the Fourier transformation. The frequency representation 
is better suited to discuss transition energies. P*^^-* has components at the sum of two frequencies of the incoming 
Ught. Since the Fourier transform of the real electric field contains positive and negative frequencies, the difference 
frequency also appears. 

If at the first step we ignore screening effects, then Eqs. (5) and (7) only contain the second-order irreducible 
susceptibility 

X^lijki^, qi, q2; t - h,h - t2) = ~ (^) Q{t - h) e{h - i2) ^ ^ ^L„/;k„r ( %l" ■,,.,AQ1z) 

xD',,^.Jq2.) [f{E^,i)-f{E^,,)] exp 



Ek„ 



--Dl||,';k||;(9i-)^k|ii";k||i'(fe) [fiE^,i')-f{Ei,,A] exp 
X exp[-rk||;";k||i'(ii -t2)] exp 



(<i - t2) exp[-rk||;';k||((ii - t2)] 

.-EkiiZ' - Ei^„l>r ' 
- [ti — 12) 



exp[-rk||;";k||!(i - ti)], 



(10) 



which is derived in App. A. We neglect the photon momenta relative to the crystal momentum. This expression, which 
forms the basis of our discussion of SFG, goes beyond the one given in Ref. [31] in that it is valid for a time-dependent 
and spatially varying laser field. Furthermore, it includes the transverse response explicitly. Equation (10) already 
exhibits the interplay between the time interval \ti — 12\ between absorptions, the photon frequencies, the dephasing 
times, and the transition frequencies. 




FIG. 4: Detailed quantum-mechanical interpretation of a process contributing to SFG. An electron is excited from a pure state 
|k||?) in tlie Fermi sea to a superposition of states |k||/) and |k||?") by the absorption of two photons. After the absorptions 
the excited electron returns to the original pure state by emission of a SFG photon at the sum frequency. The times of the 
absorptions and the emission are indicated. The heavy wavy lines denote superpositions resulting from the absorption of 
photons (indicated by thin wavy lines with arrows) at t2 and ti, while the black dots represent electrons in pure eigenstates. 



We now discuss the physics contained in Eq. (10) with the help of Fig. 4. We consider the first of the two terms 
in Eq. (10). The interpretation of the second term is similar [63]. The step functions incorporate the time ordering 
t2 < ti < t and thus guarantee causality. The system is in equilibrium until the first absorption at time t2 creates a 
superposition of the two states \\i\\l) and |k||r), denoted by the wavy line in Fig. 4. This important physics is lost in 
the interpretation illustrated by Fig. 1. The Fermi functions make sure that one of the states is initially occupied and 
the other is empty. Let us say state |k||Z) is occupied. Since the system is in a superposition of two eigenstates, the 
polarization oscillates with the frequency (i^kyi — -Ekni')/^; follows from the first exponential in the parentheses in 
Eq. (10). Such superpositions are described by the off-diagonal components of the density matrix [64]. The diagonal 
components denoting the occupation numbers of states are not changed by a single absorption. The superposition 
decays with the dephasing rate rk||j';k||( associated with this transition, making it clear why the dephasing rates rather 
than the energy relaxation rates dominate the response. A second absorption at the later time [65] ti changes the 
state into a superposition of the originally occupied state and the state |k[|Z") with its own characteristic oscillation 
frequency (ii^kii; — -EkMc)/^ ^ind dephasing rate. This oscillating polarization can emit a photon at that frequency. 
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After the emission the electron is again in the pure eigcnstatc |k||?). The nonhnear susceptibiUty in Eq. (10) contains 
a sum over many contributions of this type from different momenta and bands [51]. Note, the product of three dipole 

matrix elements appearing in XiVr is responsible for the surface sensitivity of SFG, since in inversion symmetric crystals 
the product of dipole matrix elements connecting three states vanishes except when inversion symmetry is explicitly 
broken, e.g., by the surface. 

In frequency space the nonlinear susceptibility is given by 



D 



k|J;k||i 



—Tiio + i^k||; — E]^\\i" — i?irk||i";k||; 

^'k„i";k||J'('?u)^k||i';k„ (92z) 

^ /(-5'k||i) - fi.E^\\i') 

-TujJ + ?iuj' + E^ki,;- -Ek||i' - «?irk||i';k||i 

-^k||i';k||i(ai^)^k||Z";k||Z'(fe) 



-JkO + TiUj' + El^..l' — E'kiii" — iK^is.«l":k« 



(11) 



This shows that the contribution of intermediate (virtual) states falls off with the inverse of the initial-state energy 
plus the photon energy minus the intermediate state energy, i.e., with the inverse of the detuning. This is not related 
to the lifetime broadening, but is due to Heisenberg's uncertainty principle, which allows energy conservation to be 
violated on short time scales. The frequency picture also allows to incorporate a weight factor to account for the 
frequency resolution of the detector [51]. It is of interest to note that Eq. (11) and Fig. 4 can also describe spin- 
selective electron excitations due to circularly polarized light. Including the electron spins our response theory and 
in particular Eq. (11) apply also to magnetic systems. 

To prepare the analysis of the effect of collective plasma excitations on the nonlinear optical response we now 
include the screening of the electric fields. Screening enters in two ways: First, the effective field E within the solid 
is not identical to the external field because of linear screening, which is expressed by the Fresnel formulae [25-29] 
containing the dielectric function e, which can be determined in the RPA. Secondly, the second-order polarization 
P(^) of the electron gas, which corresponds to a displacement of charge, leads to an additional electric field [66] 



El'>{r) 



)('■ 



■ - r'|5 



-r'|3 



(5(r 



r') 



Pf (r'). 



Fourier transformation leads to 



E(2)(k,i) = -47rkk-P(2)(k,t), 



(12) 



(13) 



where k is the unit vector in the direction of k. Thus only the component of P*^^' (k, t) parallel to k, i. e., its longitudinal 
part, is accompanied by an electric field E*^^^ [66], which is also longitudinal. Note that a longitudinal component of 
the electric field and of the induced polarization generally exists even for a transverse applied external field for lattice 
models [41, 42], see Eq. (11). 

Due to the linear polarizability of the solid the additional field E^^^ leads to a polarization contribution of the form 
XE^^). Since the field E^^^ in Eq. (13) is of second order in the applied field, see Eq. (5), this polarization contribution 
must be taken into account in P^^'. Doing this selfconsistently corresponds to the summation of an RPA series [31], as 
shown in App. A. Then, p/^^ oc / dti dt2 Xijk Ej Ek where now the nonlinear susceptibility x|^]. obtains an additional 
factor and is given by 



x-2(q' qi' q2; * - *i - *2) = ^ XI / 



dt 



eloL;^m(q' ^ " *) x|rr;m7fc(^' ^^^^2' t - ti,ti - tj)- 



(14) 



Here the irreducible susceptibility xU is given by Eq. (10). Since only the longitudinal component of P*^^) is accom- 
panied by an electric field, the screening factor ffjo^g appears only for the longitudinal component. This is expressed 
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by the factor q^^q^ in the explicit expression £iong;ii (q, q; t - t) = 5ij + 47rXim(q, q; * - *) Qm^j- ^iong;jm(^' q; * - is 
the inverse matrix with respect to the indices {i, qz) and (m, qz). 

This analysis is illustrated by Fig. 5. Figure 5(a) shows the diagram of the second-order susceptibility x^^^- The 
square vertex represents the additional factor of e^J^^. It is obtained from the Dyson equation in Fig. 5(b). The 

expression for Xi^r (lO) ^ called irreducible since its diagram Fig. 5(a) with the square vertex replaced by a 

normal one cannot be cut into two by severing a single photon line. 




FIG. 5: (a) Diagrammatic representation of the second-order susceptibility x'^^ in terms of the effective electric field E. P^'^' 
is the induced nonlinear polarization. The square vertex □ represents the factor e^J^^ at the frequency of P*-^'. It appears if 

one self-consistently takes into account the electric field due to the polarization P'^^ of the electron system and is given by the 

RPA series shown in Fig. (b). The wriggly line refers to the electron-electron interaction through the electromagnetic field and 
is absorbed into the matrix elements D. Note, Xiong only acts on the longitudinal field components, see text. The factor s^J^^ 
can be enhanced by the plasma resonance. 

The response theory clarifies how collective plasma excitations affect SFG. They essentially enter in two ways, both 
of which arc controlled by the full (not only longitudinal) dielectric function e: 

First, the effective electric field is expressed in terms of the external field by means of Fresnel formulae [25-29], 
which contain contributions of order 1 /e for small e. The dielectric function e becomes small if the frequency of the 
external field is close to the plasma frequency. This contribution can be interpreted as field enhancem,ent. In addition, 
the outgoing (sum-frequency) electric field Eout also contains terms that are enhanced for small e due to the Fresnel 
factors. This enhancement is most pronounced if the sum frequency is close to the plasma frequency. 

Secondly, the longitudinal component of the nonlinear polarization P*^^^ of the electron system is accompanied by 
an electric field E^^^ given by Eq. (13). Thus, the factor £j~^g appears in the nonlinear susceptibility in Eq. (14) and 

thus in P(^) [31]. This leads to an enhancement of the SFG light due to the longitudinal part of P'^^ if the sum 
frequency is close to the plasma frequency. 



B. Two-photon photoemission 

To demonstrate the similarities between SFG and 2PPE, we continue by summarizing the results of the response 
theory for 2PPE. We consider the same band structure as for SFG, which is characterized by single-electron energies 
-E-kjii- We emphasize that this band structure contains the bulk states with the z-component kz of k included in I. 

The response theory starts from the observation that the photoelectron current j{t\ k, a) of electrons of momentum 
k and spin a is given by the change of occupation of the vacuum state [kcr, out) outside of the crystal. However, in 
practice the time-dependence of ) is not measured, but only the total photoelectron yield A/'(k, cr) = J dt j{t;]i.,a). 
This is similar to SFG, where only the time-integrated intensity is measured. The response theory directly determines 
the photoelectron yield Af. To prepare the discussion it is useful to first consider ordinary single-photon photoemission. 

Single-photon photoemission: The photoelectron yield is given by 



Af(k,a)=Y^ / dtidt2%(q;ti,i2;k,CT)£'i(q,ti)£'j(-q,i2), 
q 



(15) 
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with the response function (see App. B) 



%(q;il,t2;k,cr) = ^-^7— — 2^^k<T,in;k||A(9^) 



in;k(T,ii 



X exp 



■.£^k||A -£^ka,in ■ 
^ — - — ^ {t2 - tl) 



-rk||^;k<T,in|t2— tl| 



X /(-Bk||A)£'i,A;ka,in(-'Z. 



(16) 



Here, |kcr, in) is a state with momentum k and spin a inside the crystal but above the vacuum energy. We have 
again neglected the momentum transferred by the photon. The standard diagrammatic representation of ordinary 
photocmission is shown in Fig. 6 [67]. The effective field E within the sohd should again bo expressed in terms of the 
external light field with the help of the proper boundary conditions. The response function r] will play a role when 
we discuss the various contributions to 2PPE. 




FIG. 6: Diagrammatic representation of the relation j oc J dti dt2 rjij EiEj for tlie photoelectron current in ordinary plioto- 
emission [67]. The wavy lines denote the effective electric field E within the solid and the arrow denotes the emitted electron 
current j (or the photoelectron yield), which is of second order in the electric field. The response function r] is discussed in the 
text. The dots (•) denote dipole matrix elements D. 



We briefly commend on the structure of this expression: The prefactor 7kcr,out;kCT.in/rkCT.in;ko-.in describes the prob- 
ability that electrons excited above the vacuum energy actually leave the crystal. Photocmission is often described 
by a three-step picture [68-70]: First, electrons are excited, then they are transported to the surface, and finally they 
leave the crystal. In this work we are mainly interested in the first step. The second and third steps are incorporated 
phenomenologically by effective relaxation rates Fko-aniko-an, which describe electrons dropping below E^^c before they 
reach the surface, and effective^ transition rates 7k(T,oui;k(T.in from states above i?vac within the solid to free electron 
states outside of the solid. Note, the yield is proportional to the electric field squared and thus to the intensity of the 
incoming light. 

Two-photon photoemission: The total 2PPE yield consists of the three contributions 

^2PPE ^ ^2PPE ^ ^PPE ^ ^2PPE 

corresponding to Fig. 7(a), (b), and (c), respectively. The second and third term arise from the nonlinear optical 
properties of the solid: Close to the surface the effective field E leads to a second-order polarization P*^^\ see Eq. (5), 
which is accompanied by an electric field E^^^ This field may contribute to photoemission, leading to the processes 
in Fig. 7(b) and (c). The second-order field E^^^ is also responsible for SFG accompanying 2PPE. However, this SFG 
is usually a small effect since the SFG light intensity is of sixth order in dipole matrix elements D, see Eqs. (7) and 
(10), whereas the 2PPE current is of fourth order, as is shown below in Eq. (24). This changes if the sum frequency 
is close to the plasma frequency, in which case SFG is enhanced as discussed at the end of Sec. 11 A. 

Since the diagram in Fig. 7(a) cannot be cut into two by severing a single photon line, the first term Af^f^^ is 
irreducible, while the other two are reducible. The irreducible contribution in Eq. (17) can be written as 



A/:2PPE ^ j v!J'a^{ti,t2, t3, t4) Ei{ti) Ej{t2) Ekita) Eiit^j, 



(18) 



which is of fourth order in the electric field and of second order in the incoming intensity. This is already clear from 
the simple picture in Fig. 2: The occupation of the intermediate state |2) is proportional to the light intensity. To 
reach state |3) above iJyac another absorption is required, leading to a total proportionality to the intensity squared. 
Obviously, the structure of Eq. (18) is very similar to Eq. (7) for the SFG yield: 

J(2) oc j dtdti dt2 dt3 dti x[fk{t - tl,ti - t2) xfmii " *3, ^3 " U) 

X Ejitr) Ek{t2) Ei{t3) Em{ti). (19) 
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(a) (b) 




(c) 

FIG. 7: Diagrammatic representation of contributions to the 2PPE yield jV^pp^. (a) Direct irreducible contribution A/^^f^^ 
involving four transitions induced by the effective electric field E. (b) Reducible process involving conversion of two photons 
to a SFG photon and subsequent photoemission of third order in the fields, yielding A^cda^- (c) Reducible process involving 
conversion of all four photons to two SFG photons and ordinary photoemission (of second order in the fields) induced by the 
SFG light, yielding A/'red^^- The dots (•) denote dipole matrix elements D. The square vertex (□), representing a factor of 
^longi defined in Fig. 5(b). Plasma excitations again enter through this vertex. 



Hence, we expect similar interference effects in both cases. 
The other two contributions to jV"^^^^ are 

^^^dT = -4^^ / dti dt2 dts vijl [P^^\h) Ej{t2) Ek{h) 

+ Eiih) Pf\t2) Ek{h) + Eiih) Ej{t2) Pf ^(ts)] , (20) 

and 

f^Z7 = (4^^)' / dti dt2 TH, ih) (t2), (21) 



where rj is given in Eq. (16). Since the nonlinear susceptibility x^^^ and hence P*^^) contains three dipole matrix 
elements, the reducible contributions to the photoelectron current are of higher order in dipole matrix elements 

and are thus usually small. However, the longitudinal component of P*-^-* contains a factor £i~ng. If the nonlinear 
polarization is enhanced due to a bulk plasma resonance at the sum frequency, one expects significant contributions 
from the reducible terms. The response functions rj^^^ and tf^^^ are given in App. B. 

We next consider the response functions -q, rj'^^^ , and r/^^^^ which determine the yield 7V^^^^. The bmctions t]'^''^ and 
^2PPE appearing in Eqs. (20) and (21), respectively, arc of the same general form as i], Eq. (16), but have more terms 
resulting from different orders of the time arguments. We now first present the structure of the response expression 
for the main, irreducible contribution to the 2PPE yield, Eq. (18), and then discuss its physical interpretation. Fully 
written out, Eq. (18) reads 

A/;?PP^(k,a) = ^ /"dtirft2c;t3c!i4??-/fer(qi>q2,q3;ti,t2,t3,i4;k,a)i;i(qi,ti)£;,(q2,t2) 
qiq2q3 

xEk{ci3,t3)Ei{-ci-(ii-(i2,U). (22) 
Defining the complex transition energy 

Eu.„i — Ek'j~ 

we obtain the response function 

Vijld^{<lU'i2,<l3;tut2,h,t4;K(r) = 74 ^L,m;k||Ai(9u)^k||Ai;k||A2(92z) 

fl -L kcr,in;k(7,in \ \ ^ 

Ai A2 As 



^k||J;k]|J' = ^ «rk||J;kJ|i', (23) 
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X |e(tl - is) e(i2 - ts) - U) e-'"''ll^l^''-'"(*l-*=) g-in,.||^,,^<,.„(t2-t3) g-iO..||.3;W,.n(t3-t4) 

- e(ti - t2) e{t2 - t4) e(t4 - ts) e-^'Sl'^^'^--'-^''-*'^ ^-^n^,.rM^Mt2-U) g-iOk||X,;k,|X3(f4-t3) [f{Ek^^x,)-f{Ei,,X,)] 

(24) 

There are eight terms in the curly braces, which correspond to different temporal orders of interactions with the 
electric field. Note, the dependence of 2PPE on light polarization is incorporated in the symmetries of the tensor 
Vijid^t which depend on the dipole matrix elements D. Unlike for the nonlinear optical response, these symmetries 
have not been discussed so far. It would be very interesting to determine the symmetries for surfaces of nonmagnetic 
and magnetic solids. 

|ka,in> 



|k||?ii> 



|k||?i2> 

1 1 1 — I > 

t3 t4 t2 ti t 

FIG. 8: Interpretation of one of the processes contributing to 2PPE. An electron is excited from a pure state |k||A2) in the 
Fermi sea to a pure state |kcr, in) above E^ax: by four interactions with the electric field at the times U, as expressed by Eqs. (22) 
and (24). Note, the photoelectron current is proportional to the fourth power of the electric field and thus to the intensity 
squared, as expected for two-photon photoemission. The heavy wavy lines denote superpositions of states |k||A2) and IkyAs), 
|k||A2) and |kcr, in) etc., while the black dots represent pure eigenstates. Compare with Fig. 4 for SFG. 

Equation (24) forms the basis for our discussion of 2PPE. To clarify the time dependence exhibited in Eq. (24) we 
discuss the second term, the others are in principle similar but correspond to different orders of the times ti. The 
processes described by this term are illustrated in Fig. 8. The system starts in equilibrium from the state |k||A2). 
The first interaction with the electric field takes place at time is and creates a superposition of the states |k||A2) and 
jkiiAs), leading to oscillations at the frequency (-EkyAa — E\^^^x-^)/h expressed by the third exponential factor in this 
term. The Fermi factors ensure that one of the states is initially occupied and that the other one is empty. Let us 
assume that |k||A2) is occupied. The second interaction at ti changes the state into a superposition of |k||A2) and the 
vacuum state jkcr, in), leading to oscillations at the corresponding difference frequency (second exponential factor), 
and the third interaction at t^ creates a superposition of the vacuum state and |k|| Ai). After the fourth interaction the 
electron is in a pure state above E'vac and can leave the solid with finite probability. Of course, due to the sum over 
bands there are usually several contributions of this type. Only if the superpositions decay very rapidly compared 
to the pure states, a description in terms of rate equations, as suggested by Fig. 2, is applicable [15, 39, 40]. Also 
compare the discussion of SFG above, see Fig. 4. 

While SFG is only governed by the dephasing rates but not the energy relaxation rates, 2PPE depends on both. 
This is because in the 2PPE response function t^^ppe ^j^^ change of occupation of states enters besides the polarization 
of the electron gas. whereas SFG only depends on the latter. 

Note, the 2PPE yield contains four dipole matrix elements. Thus, even for inversion- symmetric crystals parity 
does not forbid 2PPE from the bulk. However, 2PPE is sensitive to a surface region of a thickness given by the mean 
free path of electrons above -Evac- The optical penetration depth is typically significantly larger than the mean free 
path and thus does not enter here. Equations (16) and (24) also illustrate that 2PPE is sensitive to specific points 
in the Brillouin zone: The photoelectron momentum k measured by momentum-resolved 2PPE is approximately the 
same as the lattice momentum of the original unperturbed electron and also of the intermediate state due to the 
small photon momentum. These effects obviously require a theoretical description that considers the k-dependent 
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states in the solid, like our approach docs as opposed to both the random-k approximation and Bloch equations. 
In view of the importance of angle-resolved (ordinary) photoemission spectroscopy (ARPES) for, e.g., cuprate high- 
Tc superconductors, k-resolved 2PPE is expected to yield interesting results in the future. On the other hand, if 
one only measures the total number of photoelectrons, the k-space resolution is lost and 2PPE and SFG give very 
similar information. It is obvious that 2PPE has the disadvantage of being limited to frequencies wi, u)2 such that 
Ep + hoJi + fiu2 lies above the vacuum energy, unlike SFG. As stated already the SFG accompanying 2PPE is usually 
small since additional dipole matrix elements are involved. 

The response expressions show that collective plasma excitations affect 2PPE in two ways: First, exactly like for 
SFG the effective field E within the metal differs from the external field due to linear screening and is enhanced close 
to the plasmon resonance. Secondly, the reducible contributions in Eqs. (20) and (21) depend on the second-order 
polarization P(^\ which contains a factor of ei~^g, see Eq. (14). P^^^ is enhanced if the sum frequency is close to the 
plasma frequency. In 2PPE this enhancement enters only in the reducible contributions in Eqs. (20) and (21). 

III. DISCUSSION 

The aim of the present section is to discuss and illustrate the results of the response theory for time-resolved SFG 
and 2PPE. In particular, we consider time-dependent effects on the femtosecond time scale. Our results exhibit 
the intimate relation between SFG and 2PPE. We can already gain insight by studying the general structure of the 
response expressions for SFG and 2PPE. for example Eqs. (7) and (18), respectively, independently of the specific 
approximations made here. For clarity we apply our response theory to a simple model system. 

A. Time-dependent effects in SFG and 2PPE 

The response expressions of the preceding section are valid for any time dependence of the exciting laser field. The 
time enters the response expressions for both SFG and 2PPE in two ways, apart from the step functions from causality, 
cf. Eqs. (10), (16), and (24): The difference between the time arguments of electric fields appears in exponentials 
oscillating at the transition frequency of the involved electron states and in exponentials decaying with the dephasing 
rate of the superposition of the two states, and, for 2PPE, also exponentials decaying with the energy relaxation rate 
of an intermediate state. (See the discussion of Figs. 4 and 8 for the interpretation of SFG and 2PPE in terms of 
electronic excitations.) The time passing between absorptions can be controlled by the pulse shape of the exciting 
laser pulses: If the total duration T of a pulse of arbitrary shape is much larger than typical relaxation times r then 
the yield depends on the probability to absorb two photons within a time interval r, which is independent of T. On 
the other hand, for T ^ r there is almost no relaxation during the pulse. Thus the response theory reproduces the 
well-known result that r can only be inferred from SFG or 2PPE experiments if the total pulse duration is T ~ r. 

To be more specific, in most experiments two approximately Gaussian pulses are used (pump-probe method) [1-20]. 
If the two pulses are of different mean frequencies ix>\ and LO2 (two-color case) and one measures the SFG or 2PPE 
response at the sum frequency u>\ +ij02 then it is obvious which photon was absorbed out of which pulse. Then for long 
time delay AT compared to the single-pulse duration the relaxation rate of intermediate states can be read of directly 
from the AT dependence of the total yield. In pump-probe experiments with two pulses of the same mean frequency 
(jO (single-color case), photons can be absorbed out of the same or different pulses. However, the contribution with all 
absorptions out of the same pulse obviously does not depend on AT, just leading to a constant background. Note, 
in all these cases only a typical relaxation time enters, which usually is a weighted average over relaxation times of 
many states [51]. If only a single relevant intermediate state is present, e.g., for a quantum-well state, or if there 
are many but of similar relaxation rate, the relaxation time extracted from experiment will be the actual dephasing 
time of intermediate states. However, if intermediate states with very different dynamical properties are involved, 
for example if both sp and d bands are relevant, the measured relaxation time does not describe any single excited 
electron state. 

In pump-probe SHG [1-5] or pump-probe single-color 2PPE [8, 11-13, 15] experiments, time-dependent interference 
effects arc especially pronounced. Their origin is the following: The first absorption of a photon of frequency uj sets 
up an oscillating polarization of the excited electrons. Now the probability of a second absorption depends on the 
relative phase of the oscillating polarization and the second photon. Since the oscillating polarization is described 
by the off-diagonal components of the density matrix p, a description in terms of rate equations, which omits these 
components, is unable to describe interference. 

For further illustration of this interference, we show results for SHG and 2PPE for a simple model. Unless stated 
otherwise, this model consists of three bands. The lowest one is a three-dimensional tight-binding band 1 with band 
center at [71] —3.33 eV (all energies are measured relative to the Fermi energy) and half width 3.81 eV. The band 
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maximum is at k = 0. The second, rather fiat tight-binding band 2 is centered at 2.29 eV with half width 0.48 eV 
and maximum also at k = 0. Finally, there is a free electron band 3 representing electrons above the vacuum energy 
-E'vac = 4.29 eV. There exist points in the Brillouin zone for which the energy differences between bands 2 and 1 as 
well as between bands 3 and 2 both equal the photon energy of hco = 3.05 eV. We assume that the relaxation rates 
r„j„2 only depend on the band indices rii, n2 but not on the k vector (see below). We use the energy relaxation 
rates ?ir22 = 0.191 eV (corresponding to the Hfetime T2 = 3.5 fs) and HTss = 0.381 eV (ts = 1.7 fs) and no additional 
dephasing, i.e., „^ = in Eq. (4). These short lifetimes are assumed to bring out the time-dependent effects more 
clearly. The dipole matrix elements are treated as constants. 

In the following we use this model to show how time-dependent effects emerge from our response theory. For clarity 
we neglect the Fresnel formulae, which do not change the results qualitatively. We demonstrate that our theory gives 
reasonable results for a moderately complicated system. Obviously, it can be applied to a more realistic band structure 
at the expense of computation time. The boundary conditions (Fresnel factors) are also omitted for simplicity. 




5 10 15 20 25 30 
pump-probe delay AT (fs) 



FIG. 9: (a) Total yield of photoelectrons of momentum k for single-color pump-probe 2PPE as a function of the delay time 
AT betwccu pump and probe pulses. The parameters of the model are given in the text. The k vector is chosen such that 
the transition frequencies perfectly match the frequency of incoming light. Only results for AT > are shown, since the curve 
is symmetric about AT = for identical pump and probe pulses. All curves in this and the following figures are scaled such 
that the limit for large AT is unity, (b) Total SHG yield for single-color pump-probe SHG as a function of the delay time 
AT using the same parameters. The inset shows the SHG yield for flat bands with transition frequencies that match the light 
frequency perfectly. Note, for extracting relaxation rates from experimental data the resolution of the photoelectron and SHG 
light detectors must be taken into account [51]. 

In Fig. 9(a) wc show the 2PPE photoelectron yield for a particular momentum k as a function of the delay time 
AT between two identical Gaussian pump and probe pulses. A mean photon energy of hut = 3.05 eV is assumed, 
corresponding to a wave length of about A = 400 nm, and the duration of each pulse is 10.3 fs (full width at half 
maximum of the Gaussian envelope of the electric field) . The vector k is chosen so that the transition energies between 
the bands match hui. In Fig. 9(b) we show the total SHG photon yield for exactly the same system. Unlike 2PPE, 
SHG integrates over the whole Brillouin zone. Nevertheless, the overall similarity of Figs. 9(a) and (b) demonstrates 
the similarity of the response expressions for SHG and 2PPE, compare Eqs. (7) and (18), for example. It means that 
similar information, e.g., about the relaxation rates, can be obtained from both. The SHG curve is quite similar to 
the case of flat bands, shown in the inset in Fig. 9(b). This means that only a small region of k space contributes. 
The resulting interference between different k points becomes apparent in the tail of the interference pattern, where 
the main plot in Fig. 9(b) is more irregular and decays faster. This is the averaging effect discussed in Ref. [51]. More 
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precisely it is an interference effect between different oscillation frequencies of superpositions of different states. 

The 2PPE and SHG interference patterns in Fig. 9 show the well-known 8 : 1 enhancement of the signal for 
AT = 0. This enhancement is duo to the yield being of fourth order in the field: For a single pulse the signal would be 
proportional to E^, for two isolated pulses this becomes 2E'^, but for two overlapping pulses the amplitude is doubled, 
leading to {2E f = 16E^. 

For both SHG and 2PPE, the central part of the interference pattern, which corresponds to short delay times AT 
up to about the single-pulse duration T, is dominated by the four-field autocorrelation function 



j dWi CL)2 (LJ3 Ei{wi) Ej{u2) Ek{u)3)Ei{-Wi -W2- ws)- 



(25) 



This central part stems from the overlap of the two pulses and would be present even for very fast relaxation: Then 
the response functions x'^^ and t^^ppe g^j-g very sharply peaked in time and thus nearly constant in frequency space, 
leading to J^^^ oc A'-'^^ and Af^^^^ oc for the SHG and 2PPE yield, respectively, see Eqs. (7) and (18). The 
autocorrelation signal alone is shown in Fig. 10(b). 
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FIG. 10: Demonstration of the lifetime dependence of 2PPE. (a) Total 2PPE yield for the same model parameters as used in 
Fig. 9 with a lifetime of states in the intermediate band of T2 = 6.9 fs (dotted curve) and with the very small value T2 = 0.86 fs 
(heavy solid curve). The dashed curves show the exponential decay with the dephasing rate ri2 = t^^ /2 for t2 = 6.9 fs. (b) 
Four-field autocorrelation function of the pump-probe laser field. Note the similarity to the fast-relaxation result in Fig. (a). 
The black bar denotes half the laser pulse duration. 

In Sec. II we have discussed the response expressions for time-dependent SFG and 2PPE, (10) and (24), respectively. 
The first interaction creates an oscillating polarization. There is interference if the phase information is still preserved 
when the second photon is absorbed. This is governed by the dephasing time . Thus the interference effects should 
decay with the time constant for large delays AT. This is shown in Fig. 10(a) for moderately fast {t2 = 6.9 fs) 
and extremely fast {t2 = 0.86 fs) relaxation. For the slower relaxation the tail indeed decays with but to observe 
this one obviously has to look at rather large AT where the interference is already weak. For fast relaxation the curve 
is nearly indistinguishable from the autocorrelation in Fig. 10(b). 

However, there is another crucial origin of the decay of interference: Intermediate states with energies that do not 
exactly match the energy of the original state plus the photon energy lead to beats at the frequency of the detuning. 
This effect can be seen from the response expressions. We now discuss this for the case of SFG: For pulses of short 
duration T, the times ti and <2 in Eq. (10) can be approximated by the pulse centers if we are interested in phenomena 
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FIG. 11: Demonstration of the dependence of 2PPE on the detuning of the intermediate band. The heavy sohd curve shows 

the total 2PPE yield for the same parameters as used in Fig. 9 but with the intermediate band shifted downward in energy by 
0.1 hui. u! is the mean frequency of the exciting laser field. For comparison, the dotted curve shows the 2PPE yield for unshifted 
bands with matching transition frequencies. Note the beats apparent for Ar > 10 fs. This shows how detuning affects the 
decay of the signal and needs to be taJsen into account when extracting relaxation rates from experiments. 



at frequencies small compared to T ^. Then Eq. (10) shows that the polarization of the electron system shortly before 
the second interaction at ti is proportional to 

(26) 

omitting the sum over states. The last factor stems from the electric field of frequency u) describing the first interaction 
at time t2 < ti. The decaying exponential obviously describes the decay of interference with the dephasing rate 
rk||i';k||i- The oscillating terms arc of the form exjp[—i6uj{ti — t2)] exp(—iLuti) with Slo — (i5k||j' — E]s_^^i)/fi — u. Thus 
we expect slow beats with the detuning frequency Su, which lead to an initial decay of the signal on a time scale of 
{6w)~^. There should be a recurring signal at large delay times, but this is in practice suppressed by relaxation. A 
similar argument can be made for 2PPE using Eq. (24). The effect is clearly seen in Fig. 11 for 2PPE: The width of 
the pattern is reduced by the detuning. Its tails also become more irregular. 

Next, we turn to the effect of the band structure. We first discuss SHG. The SHG yield is determined by a sum 
over many transitions of diff'erent energies and dephasing rates, see Eq. (10). If the decay is governed by dephasing 
one observes the smallest dephasing rate at large time delays AT. However, at intermediate AT one sees an averaged 
rate. The dephasing rate rk||;;k||i' for two bands I and I' should usually not change dramatically with k. On the other 
hand, the contribution of detuning is necessarily different for transitions with different transition frequencies. Thus, in 
the interference pattern a continuum of beating frequencies appears. Consequently, the initial decay is governed by an 
average detuning and later, probably unobservable, recurring signals are strongly reduced by destructive interference 
of difFcront beating frequencies. The averages are weighted by a factor approximately inversely proportional to the 
detuning, as seen from Eq. (11). For narrow valence and intermediate state bands the average is restricted to a small 
effective band width W. For broader bands but constant relaxation rates throughout each band the dependence of 
numerical results (not shown) for the SHG yield on the width of the intermediate band turns out to be weak, since in 
this case all contributing processes are governed by the same relaxation rates [72]. Hence, if only a single intermediate 
state or a few states contribute significantly [48] or if there is a narrow band with uniform relaxation rates, the Bloch 
equations should work well. In this case our expressions reduce to a pcrturbativc solution of the Bloch equations. 

On the other hand, if there is strong electron-electron scattering at certain k vectors, e.g., due to Fermi surface 
nesting, the rates can be strongly k dependent. If many states of different relaxation rates enter the SHG photon 
yield, then for broad bands the description of SHG using optical Bloch equations with a single intermediate state is 
not justified. If interference patterns are fitted with results from Bloch equations, there is no simple relation between 
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the extracted relaxation rate and the dephasing rates of the excited electrons. 

For 2PPE the situation is quite different, since this method allows to probe specific momenta k in the Brillouin 
zone. Here, the band width is not crucial. There are typically several contributions to the photoelectron yield, since 
there are several unoccupied bands. The contributions are again weighted with the inverse detuning, but now there is 
only a small number of states involved for fixed photoelectron momentum k. Thus a description in terms of a small 
number of states, e.g., by optical Bloch equations, is valid. However, if one experimentally integrates over k, 2PPE 
behaves much like SFG. 

As mentioned above, 2PPE is generally accompanied by SFG, although the latter is generally smaller in intensity 
due to the appearance of additional dipolc matrix elements. It might be interesting to perform both SFG and 2PPE 
experiments on the same sample. The two techniques are complementary in that 2PPE gives information about 
specific points in the Brillouin zone whereas SFG averages over the whole zone. Furthermore, comparison of Eqs. (10) 
and (24) shows that while the general form of the expressions for SFG and 2PPE is similar, they do depend on 
the material parameters in quite different ways. We give three examples: First, 2PPE also depends on the energy 
relaxation rates (lifetimes) directly, whereas SFG only depends on the dephasing rates. Secondly, SFG crucially 
depends on the dipole matrix element D31 of the transition from the excited state above i?vac to the original state in 
the Fermi sea, whereas 2PPE does not. Thus comparison of SFG and 2PPE may proof useful for measuring the dipole 
matrix elements. Thirdly, SFG generally results from a much thinner surface region than 2PPE and the relaxation 
rates obtained from 2PPE are more bulk-like, allowing to study the dependence of the rates on the distance from the 
surface. Finally, we have shown that there is a contribution to 2PPE from SFG light generated within the solid, see 
Eqs. (20) and (21) as well as Figs. 7(b) and (c). Simultaneous measurement of SFG and 2PPE may allow to detect 
this interesting effect. 

B. Collective plasma excitations 

Since SFG and 2PPE may be strongly enhanced by collective plasma excitations, it is useful to discuss them in 
the framework of the response theory. Our goal is to show plasmon enhancement of SFG and 2PPE in principle, 
even though the plasma frequency is larger than currently accessible laser frequencies in some metals (but not, for 
example, in silver, many heavy-fermion metals, and the interesting compound MgBj). Note, the plasma frequency is 
smaller in clusters, for which the same general picture applies. 

We have seen in Sec. II that in both SFG and 2PPE field enhancement of the effective electric field E is described 
by the Fresnel formulae [25-29]. This mechanism is relevant at the frequency of the exciting laser field and, in the 
case of SFG, also at the sum frequency for the outgoing SFG light. It corresponds to the coupling of the external 
light field to plasmon-polaritons in the solid [42]. 

The second important origin of plasmon enhancement is the screening of the nonlinear polarization P'^\ which is 
caused by the effective electric field E^^^ accompanying the longitudinal part of P^^^ and appears at the sum frequency. 
Since the electric field E*^^-* is longitudinal, a true plasmon excitation is involved. It is important to remember that 
the exciting light does couple to plasmons in real solids; this coupling is only absent in simple jellium models [41, 42]. 

Note, in the case of pump-probe SFG with two pulses of diffcrc^ut mean frequencies u)\ and 0)2, one of them and the 
sum frequency Wi -|- W2 can be close to the plasma frequency. This so-called double resonance leads to a particularly 
strong enhancement [35]. Surface plasmons lie outside the scope of this paper, since they require a more detailed 
description of the surface. See Ref. [42] for a discussion. 

Motivated by 2PPE experiments on clusters [19], we briefly consider the plasmon decay. A plasmon decays into a 
single particle-hole pair [79]. The probability of this decay is determined by the phase space available for the final 
electron-hole pair. It is only energetically possible if the plasmon dispersion lies within the electron-hole continuum 
at the plasmon momentum q, leading to Landau damping. On the other hand, decay into a single electron-hole pair 
may be possible close to the surface, since translational symmetry is broken and Qz is not conserved. If the energy 
of the electron is higher than the vacuum energy, photoemission may result. Creation of several pairs is possible 
by subsequent inelastic electron-electron scattering. A plasmon also looses energy through inelastic scattering of the 
virtual electrons and holes in the loop in Fig. 5(b). This process is governed by the single-particle relaxation rates. 
The plasmon lifetime is thus shorter than typical lifetimes of the relevant excited electrons. 

A plasma mode can be multiply excited. In a recent 2PPE experiment a doubly excited plasma mode of silver 
nanoparticles decays into a single electron- hole pair [19, 80]. What is actually observed is an enhancement of the 
2PPE yield when the sum frequency is close to twice the plasma frequency. The origin of this effect is that 2PPE 
with the incident-light frequency close to the plasma frequency is enhanced due to field enhancement. The general 
process is not specific to clusters but is also relevant for flat surfaces. It would be interesting to look for this effect 
experimentally. 
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C. Further remarks 

Concerning the range of validity of the second-order response theory we remark the following. We consider first 

pump-probe SFG with a very long delay time AT. The first-order density operator p^^^ describes the result of the 
first interaction. It contains a finite polarization (off-diagonal components) but no change of occupation (diagonal 
components), see Eq. (All). A change of occupation is only obtained from p(^) and higher-order contributions, which 
involve a larger number of dipole matrix elements and are usually small compared to p^^\ However, the off-diagonal 
components usually decay faster than the diagonal ones so that for long delay times the higher-order change of 
occiipation can dominate over the second-order polarization and the second-order approximation breaks down. On 
the other hand, our expressions for 2PPE already include the change of occupation due to the first pulse, since we have 
directly calculated the photoelectron yield to fourth order. Thus the results should hold even for long delay times. 
For the case that the polarization has decayed at the time of the second pulse, but the non-equilibrium occupation 
has not, the resulting limiting form of ry^PPE jg giyg^ }-,y jrjq (B20). It only depends on the energy relaxation rates. 
This is the case where rate equations are appropriate [15, 39, 40]. Due to the vanishing polarization there are no 
interference effects. 

There is an alternative and physically appealing description of pump-probe SFG and 2PPE as a two-step process: 
The first pulse creates a non-equilibrium distribution, which is probed by the second one. We now discuss the validity 
of calculations based on this picture. In App. A we derive an expression for the linear susceptibility of an electron gas 
in an arbitrary non-equilibrium state described by the density matrix pneq, see Eq. (A18). If we insert p'^-* due to the 
first pulse for Pneq, we obtain a two-step description for x^^' and the polarization P^^^. Omitting the details, we only 
state that the result is identical to the one obtained directly for the second-order polarization P^^^. Eq. (5), but with 

the full susceptibility x^^^ replaced by its irreducible part xlrr ^Iq. (10). Thus, by assuming two separate interaction 
processes and treating each in a first-order approximation, we loose the screening of the second-order polarization. 
This is not justified if the sum frequency lies close to the plasma resonance. 

Next we consider a two-step description of 2PPE: The 2PPE photoelectron yield ^2PPE expressed in terms of 
an arbitrary non-equilibrium density matrix pneq as discussed in App. B. Then the second-order density matrix p^^^ 
due to the first pulse is inserted for pncq- We reobtain the full irreducible fourth-order result Af^^^^ of Fig. 7(a), but 
only part of the reducible contributions. Fig. 7(b), (c): The two-step description neglects contributions of two photons 
out of different pulses being converted into one SHG photon. These contributions may become important if the sum 
frequency is close to a plasma resonance. In conclusion, the two-step picture of SFG and 2PPE is valid unless the 
response at the sum frequency is enhanced by plasmon effects. 

Finally, we emphasize that our theory can also describe timc-rcsolved SFG and 2PPE from ferromagnetically ordered 
systems. Ultimately, the light couples to the (spin) magnetization through spin-orbit coupling, which is incorporated, 
in principle, in the dipole matrix elements D and the band structure. The spin-dependent matrix elements can be 
calculated by a perturbative expansion in the spin-orbit coupling [31, 47]. SFG and 2PPE also depend on magnetic 
order through the band energies E'kyi and relaxation rates rk||;;kj|;' > since / also contains a spin index a. Of particular 
importance for magnetically-ordered materials is the rotation of the polarization of SHG light relative to incident 
light (NOLIMOKE) [31, 47, 48]. As mentioned above, the light polarization is controlled by the symmetries of the 

tensor Xi/ki which are known for low-index surfaces [78]. 

Compared to NOLIMOKE, 2PPE for magnetic systems has the advantage that in principle one can obtain infor- 
mation on the spin-dependent lifetimes of electrons in specific states |k||Z). The dependence of 2PPE on the light 
polarization for magnetic systems has not been studied so far. In the response theory this dependence is controlled 
by the symmetries of the tensor rjfji^i^, as mentioned in Sec. II B. 

For pump-probe experiments with long time delay AT the main contribution to 2PPE comes from the change 
of occupation brought about by the pump pulse. Only in this case the photoelectron yield is proportional to the 
occupation of the corresponding intermediate states. If in addition the matrix elements and the relaxation rates out 
of vacuum states depend only weakly on spin, then Eq. (B8) shows that the 2PPE yield becomes proportional to the 
spin-dependent occupation of these intermediate states: 

A/-2PPE(Ar; k, <j) a Pncq;k.a;k..(AT) = nncq;k,.a ( AT) , (27) 

where n^ccx-Mva denotes the non-equilibrium occupation of the state [kz^a) after the pump pulse. Then the difference 
of the spin-up and spin-down 2PPE yield, A/'2PP'^(AT; k, t) - A/'2PPJ^(AT; k, [), is proportional to the difference of 
the occupations and thus to the transient magnetization of the intermediate states. 

Note, circularly polarized light might excite electrons spin-selectively due to angular-momentum conservation. In 
our response theory these selection rules are incorporated in the dipole matrix elements D. Conversely, spin-selective 
excitation will lead to corresponding polarization of the SFG light. The use of circularly polarized light in 2PPE and 
SFG is of particular interest regarding ferromagnets and transient magnetizations. 
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D. Conclusions 



To summarize, we have presented a unified perturbative response theory for time-resolved SFG and 2PPE. The 
theory is fully quantum-mechanical and contains the interference effects described by off-diagonal components of the 
density matrix. It does not rely on any assumption about the time or frequency dependence of the exciting laser 
pulses. The solid is described by its band structure, relaxation rates, and dipole matrix elements. We have discussed 
metals but the response theory can be applied to semiconductors and insulators as well, see, e.g., Ref. [81]. Since 
the theory is formulated directly in the time domain, it presents a suitable framework for the discussion of the time- 
dependent physics of SFG and 2PPE. We have shown that similar information as from 2PPE can be gained from 
SHG. Of course, 2PPE is sensitive to specific momenta k in the Brillouin zone, while SHG in general is not. A simple 
tight-binding model of a metal has been studied in order to show that the theory gives reasonable numerical results 
and to illustrate the following effects important for the understanding of SFG and 2PPE. 

We have shown how relaxation rates and detuning affect the interference patterns in single-color pump-probe SHG 
and 2PPE experiments: The lifetime in the intermediate states and their detuning with respect to the photon energy 
lead to a similar narrowing of the interference patterns. The effect of detuning must be taken into account in order 
to extract meaningful life times from such experiments. Also, in particular in SHG the measured relaxation rate 
is a weighted average over the relaxation rates of many excited states. Furthermore, the weights in this average 
change with the pump-probe delay. Thus different rates govern the decay of the interference pattern depending on the 
pump-probe delay — the decay is not simply exponential. We have also discussed the range of validity of the optical 
Bloch equations, applicable if only a few states contribute, and of semiclassical rate equations valid for very long 
pump-probe delays. Both approaches are limiting cases of our theory. 

Finally, we have considered the role played by collective plasma excitations. Plasmon effects in both SFG and 
2PPE can only partly be understood in terms of field enhancement at the surface. One also has to take the electric 
field accompanying a nonlinear polarization of the electron system into account. This effect leads to interesting 
additional contributions to 2PPE, in which incoming photons are converted into sum-frequency photons which then 
lead to ordinary photoemission. These contributions should be observable if the sum frequency is close to the plasma 
frequency. 
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APPENDIX A: RESPONSE THEORY FOR THE NONLINEAR OPTICAL RESPONSE 

In this appendix we derive the transverse second-order susceptibility and polarization for arbitrary pulse shapes 
of the exciting laser field. The resulting expressions allow to calculate the SFG yield for arbitrary pulse shapes, 
thereby going beyond the results of Hiibner and Bennemann [31] for continuous-wave, monochromatic light. The 
self-consistent-field approach [30] is applied to a solid described by its band structure and relaxation rates. The flat 
surface is assumed to lie at ^; = with the solid at 2; < 0. We neglect the intraband contribution, which is reasonable 
for optical frequencies. 

The single-particle Hamiltonian is H = Hq + V, where Hq describes the unperturbed solid with the normalized 
eigenstates |k||/) and eigenenergies Skyj. ky is the crystal momentum parallel to the surface and all other quantum 
numbers, discrete as well as continuous, are collectively denoted by I (see the discussion at the beginning of Sec. II A). 

The time-dependent perturbation is [74] 

with the vector potential A, which is treated classically. We have made the usual approximation to neglect the 
quadratic term in A and have used a gauge with vanishing scalar potential. We will later need the temporal Fourier 
transform (using the convention of Ref. [31]) 



V{r,w) = J *l/(r,i) 



— E(r,u;)- V+-^[V-E(r,u;)], (A2) 
mu 2mu) 
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where we have used E = — (1/c) dA/dt. Inserting the spatial Fourier transform we get 



^(r' ^^ = -^12 • ( V - y ) • (A3) 

q 



The matrix elements of V are 



(k||/|y|k||+q||,0 = -^^E(q,a;)-(k||^|e-^'i-"-(V-zq/2)|k||+q||,0 



e X] ^(^' ^) ■ i;k|| +q|| (g^ , w) (A4) 



9z 

using momentum conservation and q = {ci\\,qz). Here, 

Dk„i;k„+q,„i'fe,u;) = ^ (ky^le-^'i-"- (V - zq/2) |k|| + q||,r). (A5) 

If the field E were purely transverse the term jq/2 in the parentheses would drop out, but this is not guaranteed close 
to a surface. It is not our goal to calculate D explicitly. We only remark that if one uses the dipole approximation 
[26-28, 47 50, 57, 73] e^^'^ = 1, the contribution from iq/2 vanishes since we neglect the intraband contributions so 
that V ^l, and the remainder gives [62, 75] 

A (k||ZlV|k,i +q||,0 = -i- (k||?|[i/,r]|k„ +q„,/') 

= - ^'"'"^1'"^""" ('-II^Mkll+qihO- (A6) 

Since the response is dominated by contributions with Tiuj ~ E^ky+qn ,i' ~^is.\\l the prefactor can be further approximated 
by unity if the frequency spectrum of the incoming light is sufficiently narrow. The dipole approximation should be 
justified since the electric field changes slowly on the scale of the lattice constant (the skin depth is about one order 
of magnitude larger than the lattice constant). However, our response theory for SFG does not require the dipole 
approximation to be made. 

The time evolution of the density operator p is described by the master or von Neumann equation [43, 61] 

The functional 7?.[p] represents relaxation terms made explicit below. Matrix elements of p are written as pky i;k]| = 
(k||Z|p|kj|Z'). The master equation then reads [43] 

^Pk||i;kj|i' = ^(k||i|[J?o + l^,P]|kj|^') + ^k||kl|<5H' 

X / , 7k||i;k]|'i"Pk]|';";kj|'i" - rk||i;kj|i'Pk||i;kj|i'- (A8) 

Here, rk||;;k||i = t^"; is the inverse lifetime of state |k||Z), which arises mainly from inelastic electron-electron scattering. 
7k||;;k]|;' gives the rate of spontaneous transitions from state |kj|/') to state |k||/). Because of conservation of electron 
number 

^k||i;k||i = 2^ 7kj|i';k||;- (A9) 

k||Z' 

Primed sums run over all states except |k||Z). rk||i;k||( and 7k||(:k]|(' describe energy relaxation, i.e., the change of 
the diagonal components of p, whereas the dephasing rate rk||;;k]|(' with |k||^) ^ Iky/') describes relaxation of the 
off-diagonal components. 
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To solve the master equation (A8) perturbatively, the density operator is expanded in powers of the perturbation 
V as p = p^"-' + p*^^^ + p'-^^ + . . . In thermal equilibrium the unperturbed density matrix is expressed in terms of the 
Fermi function, p^°Vk' i' ~ ^kyk' Sw f{Ek^^i). The temporal Fourier transform of Eq. (A8) reads 



«Wpk||i;k'i'(w) 



^^k||i -^kj,;' 



rk||J;k||i'J Pk||(;k||;'(w) 

^kl ^e-^"'E [(k||^ln^)lM^OPk;,'r;k|,H^)-Pk,,^;k|,'r(^)(k^[^1yWI^0 



+ "^kiikjif^H' X! 7k||;;k]|';"Pk||'i";k||'('A3^) 
k"i" 



Keeping only terms linear in V one obtains 

Pk||U||+q||,i'H = ^ Z]°k||(;k||+q||,i'(fe>'^) " E(q, w) 

/(i^k||+q||,iO-/(igk||0 
— TiW + i?k||+q||,i' — -Ekiii + «?irk||;;k||+q||,;' 



(All) 



Note that the diagonal components vanish: There is no change of occupation to first order. 

The polarization is given by the thermal average of — eD, which is the conjugate of the electric field according to 
Eq. (A4), 



P(q,w) = -^^^Pk||;;k||+q||,;'MDk||+q||,(';k||;(-92>'^) 



where v is the volume. To first order 

P;^^^ (q|| , , w) = ^ (q|| , 9^ , 9^ , w) (q|| , , w) , 

I'z 

with the linear susceptibility of Lindhard form 

2 



(A12) 



(A13) 



k|| W 



X'DL;;k,|+q||,('(9z,'^) 



/(■E'k||+q||,i') - /(i^k|| 



—?ko + i^kii+q,,,;' — -Bk||Z + i?irk||Z;k||+q||,i 



(A14) 



shown diagrammatically in Fig. 3. It takes into account that the z component of momentum is not conserved. 

If we neglect the frequency dependence of D we obtain a result in the time domain. Equation (A13) and the Fourier 
transform of Eq. (A14) give 



xEj{q\\,q'z,ti), 



(A15) 



with 



27ri 



(q||, 9z, Qz,t-tl) = -^Q{t-h)Y,Yl ^k||+q||,i';k||i(-9z) D 



X [f{Ek +q..,i')-f{Ek i)] exp 



k||i;k||+q|| 



k|| W 

•-^k||+q||, 



n 



-it-h) 



exp[-rk||i;k||+q||,i'(i-ii)]- 



(A16) 



The response is thus only non-zero if i > ti , which expresses causality. 

The above results have been obtained under the assumption that the system is initially in thermal equilibrium. 
We now consider the non-equilibrium case. Then the unperturbed polarization P^eq is, in general, non-vanishing so 
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that the electrons experience an effective field even in the absence of an external perturbation, leading to a master 
equation that is nonlinear in the non-equilibrium density operator p^eq- We assume, however, that this effect of 
electron-electron interaction is negligible. Then the linear response of a non-equilibrium system can be written as 



/oo 
d^' Xne<i;ij (^11 ' 9^' '^5 ^11 ' 9^, w') Ej (q' 
-oo 



a;'), 



(A17) 



with 



xiJq;ij{<l\\,qz,u};q;pqi,u)') = -- ^ ^k||+q||,i';k||i(-fe,'^) 



^k||i;k||+q' ;''(9i.'^0pi°eU||+q;,i'^k||+q||,i'(^ - 



neq;ki|;;ki|+q|i-q; 



, ,„{ij-uj')Dl 



+q||-q' i";k||+q||,; 



(A18) 



This equation gives the linear susceptibility in terms of an arbitrary density operator pneq- 

To return to the response of an equilibrium system, we now consider the second-order contribution. We collect the 
terms in the master equation (AlO) that are of second order in the effective electric field. 



■ (2) 

*'^Pk||J;k||+q||,;' 



1 

ih 



k||i -^k||+q|| J';Pk||/;k||+q||,i' 



(2) 



+ 1 (k|| 1 1 [1^(1) , |k|| + , r) + 1 (k|| / 1 , |k|| + , r 



ih 



r r (2) Ti (2) 

+ 0q||00H'2^ 7k||!;K||APK||A;K||A -ik||;;k||+q||,i'Pk||i; 



k||+q||,;'- 



(A19) 
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y^^^ =V\s the perturbation by the effective field. Higher-order perturbations n > 2, result from the electric 

field due to the displaced charge calculated at order n. From the expression for the electric field due to a polarization 
P("), Eq. (12) [66], one obtains E(")(q) = -47rqq • P("Hq)> where q = q/|q|. Importantly, this additional field also 
has to be taken into account as a perturbation [31]. Specifically, the longitudinal part of P^"' leads to a perturbation 
with matrix elements 



(k|| |k|| + q|| ,l')=eJ2 E^"^ (q, u;) ■ Dk„;;k„+q„ ,i' {q.,io), 



(A20) 



see Eq. (A4). Note, due to the reduced symmetry at the surface a transverse electric field in general leads to a 
polarization with a longitudinal component. 

Since the field E^^) = — 47rqq • P^^^ is explicitly of second order in the applied field, see Eq. (5), it must be taken 
into account in our calculation of the second-order response. Prom Eqs. (A12) and (A19) we then obtain for qii ^ 



k||+q| 



.P:k„;(-9^''^) 



X / d\'Y.j rf-'[^^,,,k,,+q,,,A(«i>-OpL;Vq;,,A;k,,+q,,,.(---') 



A^k||i;k||+q||-q' A 



k||+q||-q;,A;k||+q| 



+ 477 [/(Sk„+q,„;0 - /(i^k„0] / dg;i5(||,.k„+q„,,,(g;,a;)g,-gfePf (q',a;), (A21) 

where in the last term qj| = q||. Note that the nonlinear polarization P^^^ appears on both sides of the equation. 
Solving this equation for P*^^) we obtain 



pPM=Y, / d^'xgi(q,q';a;,a;')£;,(q',a;')i;fc(q-q',a;-a;'), 



(A22) 
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with the second-order susceptibihty 

X-2(q,q';'^,'^') = ZlZl ^loig;im(q|h9^'«^.^)xS„,•fe((q||,Kz),q^q-q';w,a;') (A23) 

and 

i \ X , A ( \ (qii''^^)fe (qih«^^)j /Ao^^ 

Kqih'^zJI Kqih'^zJI 

Here, ^i^^-ij (q|| , 9z, ^2,0;) is the inverse matrix of £iong with respect to the indices {i, q^) and [j, k^). Equation (A24) 
means that eiong only acts on the longitudinal component and is unity for the transverse ones. Solving for P*-^-* 
is found to be equivalent to the summation of an RPA series for the electron-electron interaction mediated by the 
electromagnetic field. In this language the interaction is absorbed into Xim in Eq. (A24). 
The irreducible susceptibility in Eq. (A23) reads 

.(2) ^„ „/ /A_ e^Y^Y^ ^k||+q||,;;k||;"(-9^:'-') 



xLr;i,fc(q'q'>q";'^>'^') = -—J2J2 



^ Wl/' ■^'^11+111'' ~ ■^'^11'" + *^rk||Z";k||+q||,; 

'j M Mnfe (" '^ /(^k||+q||,0-/(igk,|+q;,.i') 

^^k,r;k„+q;,,.('/.,- ) i^k„+q;,..;k„+q„,fe,- - U.)_j^^ ^, ^ ^^^^^^^^^^ _ ^^^^^^^^^^ ^ Z;.rk„+q;,,.;k„+q„, 

, , k II I /(^k||+q||-qj|,;') ^ ./(-EkiiZ") 

- ^'k„ +q,, -q;,.;k,, +q,, > ) ^k,, Z^k,, +q,, -q;, . - ) r^j;^^ 

shown diagrammatically in Fig. 5(a). Finally, the time dependence of the polarization P^^) is obtained by Fourier 
transformation of Eq. (A22) using Eq. (9), leading to Eq. (10). 

APPENDIX B: RESPONSE THEORY FOR PHOTOEMISSION 

In this appendix we give details of the derivation of the time-integrated photoelectron yield. We also present the 
analytical expressions for the response functions omitted in Sec. II B. The starting point is again the master equation 
(A7). The terms of order n > 1 can be calculated recursively, 

_^ (n) -V \ 

^^^k||i;k||+q||,i' - ^^'^^•'ll' -^k||+q||,i'jPk||;;k||+q||,i' 

1 " 

+ ^E(k||^l[^^'"^p^"-'"^]|k||+q|hO 

II A;K II A 

- rk||i;k||+q||,;'Pk||i;k||+q||,i'- (^1) 

Here, y^™) is the perturbation of order m, see the discussion leading to Eq. (A20). Among the states Ik||i) etc. 
appearing in Eq. (Bl) are states |kcr, in) lying in the crystal above the vacuum level. We assume that electrons leaving 
the crystal are in states |k(T, out) and are detected without further interaction and without returning to the solid. 
Then the only way their occupation can change is through spontaneous transitions out of jkcr, in), governed by the 
rate 7k(T,out;k<7,in- Note, in principle higher vacuum bands appear by shifting the (nearly) free electron dispersion back 
into the first Brillouin zone. We omit these bands for notational simplicity. 

First, we consider the irreducible part /o[^\ This is the contribution of only the direct, first-order perturbation 
V'^^^ = y at every step of the recursion. The resulting equation reads 

^Pirr;k||i;k||+q|| ,;'(*) = ^^"^"^^ll' ~ "^•'ll +111 ) ^irriky Z;k|| +q|| ,i' (*) 

+ E ;;k|| +q;| , A (gD /'itikiiVqli , A;k|| +q|| ,i' (*) " /^telkiiUn +q|| -qji ,A (*) ^k,, +q|| -q; ,A;k|| +q|| ,V (^^ )_ ' E(q', t) 

q'A 

+ ^qiiO^Zi'Xl Tk||(;K;||ApL"U||A;K||A(*) " rk||i;k||+q|| pL"!k||;;k||+q|| ,i' (*)• (B2) 

K.|l A 
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We assume D to be frequency-independent, see App. A. In Eq. (B2) we write the first sum in the form X^q'A ^ 
reminder that all components of the external momentum q' are summed over. Hence, we here exclude q'^ from A. 
The solution for the off-diagonal elements is 

ft 



(n) 

/'irrikiiiikii+qi 



q'A 



e W^'^w^-^w exp 



k||+q||./' 



h 



{t - h) 



Dk||;;k||+q;|,A(9i)p£;ki|Vq;i,A;k„+q||,Z'(*l) -pL"r;4^^^ • E(q', *l)-(B3) 

For the diagonal components there is an additional contribution from the relaxation of secondary electrons into the 
given state out of higher-energy states [39, 40], i.e., the term with qy = 0, 1' = I in Eq. (B2). The diagonal components 
can be written in the implicit form 



In) 
Arr;k| 



^irr;k||-Kq' A;k||iV^lJ 



'=>^/*) = ^ /* Ee"'"'"""'^*"*^^ [Dk„,k„+q,,,A(9Dp[; 

q'A 

-Ptv;k,,-q' A(*l)Dk,,-q|,,A;k,,i(9;)] •E(q',il)+ / ciilE'^k„i;K„ApL".U„A;K„A(*l)- (B4) 

If the contribution of secondary electrons is small, Eq. (B3) also applies for qy = 0, I' = I. 

The photoelectron yield for momentum k and spin a is given by the time integrated photoelectron current. Equiv- 
alently, it can be written as the occupation of the appropriate vacuum state, 

A/'(k, a) = Pko-,out;k<T,out(i oo), (B5) 

since electrons are assumed not to leave the states |kcr, out) again. For the electron states outside of the crystal 



Pk(T,out;kcr,out — 'Tkcr,out;k(T,in/^k(T,iii;kcr,in) 



(B6) 



since their occupation can only change due to electrons leaving the crystal. With Eq. (B5) we find 

/oo 
dt Pi^a,in;ka,in{t). (B7) 
-oo 

The irreducible contribution to order n is obtained by inserting the irreducible part of p^"'*, Eq. (B4), into this 
equation. After changing the order of integrals we obtain 



/I \ Tk(T,out:k(T,in / j± 

^^'"^^ih Fw • k ' / ^ 

k(T,in;kcr,in J —oo 
Dk(T,in;k||-|-q||,A((i'z) Pirr;k|| +q|| , A;k(T,in (^) 



irr;kcr,in;k|| — q[ 



.aW Dk|| -qn ,\--kaMi<l^-) ■ E(q, t), 



(B8) 



where we have used that there is no relaxation into the states above i?vac in the solid. Note, Eq. (B8) describes 
photoemission out of any (possibly non-equilibrium) state. 

For ordinary photoemission, Af^^\ we have to calculate the density operator to first order, p^^K which is purely 
irreducible. Consequently, the full ordinary photoelectron yield is obtained by inserting Eq. (B3) for n = 1 into 
Eq. (B8), 



A/'(2)(k,CT) =E j dhdt2^l^J{c^■,tut2■,k,a)Ei{cl,tl)Ej{-c^,t2) 



(B9) 



with the response function 
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It is useful to write our results in the frequency domain. For ordinary photoemission this yields 

/oo 
doj r],j (q, w; k, a) E,{q, w) Ej{-q, -oj) 
-OO 



(Bll) 



with 
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The third-order contribution, Af^^\ is of interest since the irreducible third-order response appears in the reducible 

(3) (2) 

contributions to 2PPE. To calculate A/j^r from Eq. (B8), we need the off-diagonal elements of pl^^ only, i.e., the 
polarization of the electron system, which can be obtained from Eq. (B3) alone. The result is 



Xrr^(k,(^)=^ J dtldt2dt3vSli^,<l';tl,t2,t3;-k,C7) 
qq' 
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Here, fiki;k'i' = (^^kz — E\^iii)/h — iT\a;V.'i' is a complex transition frequency. The four terms in Eq. (B14) correspond 
to different time orders of interactions with the electric field. In itself, TV^^) is usually negligible compared to TV^^^ 
for the following reason: The three frequencies of incoming photons have to add up to zero so that one has to be the 
negative of the sum of the other two. However, then the sum frequency is already present in the exciting laser pulse 
and ordinary photoemission dominates the signal. 

Finally, the irreducible contribution to fourth order has the general form 



E / 



dt\ dt2 dts dti 



X vlWi^^ q'> q"; tiM,h, U; k, a) Ei{q, tx) Ej (q', is) 

x£fe(q",t3)£;i(-q-q'-q",i4). 
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jj2PPE pg^jj found by inserting Eqs. (B3) and (B4) into Eq. (B8). The photoelectron yield is determined by the 

off-diagonal components of p^^\ which in turn depend on all components of p^'^\ including the diagonal ones. New 
physics enters here: The 2PPE current depends on both the polarization and the change of occupation to second 
order. 
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If the increase of the occupation due to secondary electrons is small to second order, we can use Eq. (B3) to calculate 
all components of p'^^-' . The change of occupation due to dipole transitions and to relaxation out of excited states is 
included in Eq. (B3). Then the response function ry^^^^ reads 

4 

2PPE / f ff 4. 4. 4. 4. \ \ ^ Tk(T,out;kcr,in \ ^ j^i / \ j-^j / / \ 

Vijkl (,q>q>q ;*l,*2,*3,l4;K,Crj - -j 2^ ^ka,in;k||+q||,AW«;^k||+q||,A;k||+q||+q' A'Wz) 

il k(T,in;kcr,in yi " 
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- - U, U - t2, t2-h; 1, 0; 1, 3; 2, 3) + F{ti - U, ti - t3,t3 - t2; 1, 0; 1, 3; 1, 2) 

- F{t4 -tuh- t2, t2-t3; 0, 3; 1, 3; 2, 3) + F{ti - ti, ii - ta, - ^2; 0, 3; 1, 3; 1, 2) 

+ F{t4 - is, t3-ti,h-t2; 0, 3; 0, 2; 1, 2) - F{U - ta, is - t2, t2 - ii; 0, 3; 0, 2; 0, 1)] (B16) 
with the auxilliary function 

F{Ati,At2, Ats; ni, riz; ns, nr, ng, ng) = e(Aii) e{At2) Q{At3) 

X e-inni.n^Atl g-ian3.n4At2 g-iOn^.n^Ata [/(E^J - /(£„J] , (B17) 

where the states \ni) are defined as 

|0} = |ka;in}, |1) = |k|| + q,, , A), 

|2) = |k||+q||+qj|,A'), |3) = |k||+q||+qj|+q|[,A"). (B18) 

Thus the first term in the brackets in Eq. (B16) reads 

e{t - h) e{h - t2) Q{t2 - is) e-^"''ii+-'r^^''-'"(*-*^^ 

^e-'''^^^-'^^^-'^^-'^^^^^^^^^^^^^^ (B19) 

=0 

etc. rf^^^ has the same structure as only with more terms due to more possible time orders. One should exclude 
terms from Eq. (B16) that correspond to processes for which the system returns to the equilibrium state after two of 
the four interactions. These processes only contribute a small correction to the numerical prefactor of the ordinary 

photoelectron yield. 

If, in addition, the typical time scale of the experiment, e.g., the delay time in the pump-probe case, is long 
compared to the dephasing times, 2PPE can be described by the change of occupation alone. Interference effects 
are then absent. In this limit, only the terms F in Eq. (B16) with = 714 contribute, where ns is an excited state 
reachable by a single interaction out of the Fermi sea. Then jy^PPE gjixiplifies to 

4 

2PPE/ \ ^ Tko-jOutikcTjin V~~^ 

WjfeT (q> q > q ; *i> *2, is, tr, k, cr) = — ^ — '—^ ^ 

Tl J- k(T,in;k(T,in 
^ '^qjl +q;|' ,0 -Okcr,in;k|| +q|| , A (<?^ ) ^ky +q|| , A;k|| +q|| +qj| , A' (^z) 
X -Ok||+q||+qJ|,A';k||+q||,A(~9i) -Dk||+q|| ,A;ko-,in(-9«) 

- F(ii - i4, i4 - i2, i2 - is; 1, 0; 1, 1; 2, 1) 

+ F(ii - i4, i4 - is, is - i2; 1, 0; 1, 1; 1, 2) 
- F(i4 - il, il - i2, i2 - ia; 0, 1; 1, 1; 2, 1) 

+ F(i4 - il, il - is, is - i2; 0, 1; 1, 1; 1, 2)] . (B20) 

Note, T^^PPE j^Q^ proportional to exp(— FnAi) from the second exponential in Eq. (B17), where Fn = Tj"^ is the 
energy relaxation rate of state |1) = |k|| + qy , A) and At is the time between the second and third interaction. This is 
easy to understand: After the second interaction the electron is in the pure state which decays with the rate Fn. 

Reducible contributions to the photoelectron current result from nonlinear optical effects in the solid. They are 
obtained by replacing the effective electric field E in Eqs. (Bll) and (B13) by the electric field to second order, see 
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Eq. (13). For 2PPE wc obtain the contributions given in Eqs. (20) and (21) and shown diagramniaticahy in Fig. 7(b) 
and (c). There are also contributions given by ry times a product of E and the third-order polarization P'^^. These 
are only significant if the incident light contains a frequency component large enough to allow ordinary photoemission 
and we neglect them. 

Finally, it is also possible to describe (ordinary) photoemission out of a general non-equilibrium state described by 
the density matrix pneq- The irreducible contribution is obtained from Eq. (B8) for the photoelectron yield Af^"^ by 
expressing /o|^~^'' in terms of the lower-order pl^^'^^ with the help of Eq. (B3). This is exact, since only off-diagonal 
components of p[^~^^ are needed. Then p^""^^ is replaced by Pneq- There is also a reducible part: Two photons can 
be converted into a single one at the sum frequency, which for non-equilibrium, when a finite polarization exists, can 
lead to a change of occupation of states above Svac and thus to photoemission. 
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